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1 Introduction 



It would be difficult to find a more ubiquitous polynomial family than the 
Hermite polynomials H n (x), defined by the Rodrigues formula 

H n (x) = (-1)" exp (x 2 ) — exp (-x 2 ) , n = 0, 1, 2, . . . . 

They appear in several problems of mathematical physics the most 
important probably being the solution of the Schrodinger equation [Ti] . 
Being the limiting case of several families of classical orthogonal polynomials 
[To] , they are of fundamental importance in asymptotic analysis |33| . 
The Hermite polynomials satisfy the orthogonality condition 

oo 

J e~ x2 H m (x)H n (x)dx = ^fu2 n n\b mn . 

— oo 

the differential-difference equation 

H n+1 + H' n = 2xH n , (1) 

and the reflection formula 

H n (-x) = (-l) n H n (x). (2) 

The zeros of the Hermite polynomials have several applications, notably 
in Gauss' quadrature formula for numerical integration [T^], [2H], jHO]- Several 
properties and their asymptotic behavior were studied in [T], [2], JZj ; 
Hi EE] EH and [32]. 

The asymptotic behavior of H n (x) was studied by M. Plancherel and 
W. Rotach in |2Z| using the method that now bears their name. F. W. J. 
Olver obtained asymptotic expansions for the Hermite polynomials as a 
consequence of his WKB analysis of the differential equation satisfied by the 
Parabolic Cylinder function D u (z) , related to H n (x) by 

H n (x) = 2t exp f^j D n (V2xj . 

A similar analysis using perturbation techniques was carried on by A. Voznyuk 
in 
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As an application of the results from his doctoral thesis on the multiplication- 
interpolation method, L. Heflinger CHj established asymptotic series for the 
Hermite polynomials. In [39], M. Wyman derived asymptotic formulas for 
H n (x) based on one of their integral representations. 

In this paper we shall take a different approach and analyze the differential- 
difference equations that the Hermite polynomials satisfy ([TJ using the tech- 
niques presented in ^Uj. A similar method (which we may call the discrete 
WKB method) has been applied to the solution of difference equations [S], 
jH], [12], |HE] and it is currently being extended [TT] . |3~5j . j3*o] . |3*7j . to include 
difference equations with turning points. Another type of analysis, based on 
perturbation techniques, was considered by C. Lange and R. Miura in [T7j . 
CHI, d, U EH, and m 

2 Asymptotic analysis 

We consider the approximation 

H n (x) ~ exp [f(x, n) + g(x, n)] , n — ■> oo (3) 

where 

9 = o(f), n -> oo. (4) 
Note that since Hq(x) = 1, we must have 



and 

Using (jSJ) in (JTJ), we have 



f(x,Q)=0 (5) 
g{x,0) = 0. (6) 



exp ( f + ^l + l^l + g + ^l] ( 7 ) 
\ dn 2 dri 2 dn J 



where we have used 



df 1 d 2 f 

f(x, n + 1) = f(x, n) + — (x, n) + - (a?, ra) + 
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Simplifying (J7J) and taking (jlj) into account we obtain, to leading order, 



, ld 2 f dg\ dg ( df 
exp — — — H — — J H — — exp - 



Using JHJ) in (0) we get 

_ 2 dn 2 dn J ' dx 1 V 

or, to leading order, 

1 d 2 f dg dg ( df 



2 dn 2 ^ dn ^ dx ^ \ dn 
2.1 The ray expansion 

To solve (jSJ) we use the method of characteristics, which we briefly review. 
Given the first order partial differential equation 

F(x,n,f,p,q) = 0, 

where 

df df 
dx dn 

we search for a solution f(x, n) by solving the system of "characteristic 
equations" 

dx dF dn dF 

dt dp ' dt dq' 

dp _ dF dF dq _ dF dF 

dt dx ^ ' df " dt dn ^ df 1 

aJ_dF_ + dF_ 

dt ^ dp ^ dq' 

where we now consider {x, n, f, p, q} to all be functions of the variables t and 
s. 

For (jHJ), we have 

F (x, n, f, p, q) = e q + p — 2x (10) 
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and therefore the characteristic equations are 



dq 

2 — = 
dt dt dt dt 



and 



dx dn dp 

df 



dt 



p + qe 1 



Solving (JTTJl subject to the initial conditions 

x{0,s) = s, n(0,s) = 0, q{0,s) = A(s), 

we obtain 

x = t + s, n = te A , p = 2t + 2s-e A , q = A, 
where we have used 

= F| t=0 = e A +p(0,s)-2s. 



From (JHJ) and (|TH|) we have 
which implies 

d 



/(0,s) = 0, 



= -^/(0,s) 
ds 



df dx df dn 
dx ds dn ds 



t=o 

A 



= p(0, s) x 1 + q(0, s) x = 2s - e y 
Thus, A(s) = ln(2s) and (jUj) becomes 

x = t + s, n = 2ts, p = 2t, q = ln(2s) 



(11) 
(12) 

(13) 

(14) 



(15) 



(16) 



with t > and s > 0. Since s > 0, we shall consider only the region x > 
for now. Using (|16|) in (|12|) and taking (|15|) into account, we obtain 



f(t,s) = t 2 + 2sln(2s)t. 

Solving for t and s in terms of x and n in (jlfi)) . we get 

x 1 x 1 

t — — ± -<J, s = — =F _ cr 
2 2 ' 2^2 



(17) 



(18) 



with 



a = V x 2 — 2n. 



(19) 



For a to be a real number, we shall impose the condition x > \J1n. Since 
(for a fixed value of n) we have t — > as x —>■ oo, we consider the solution 



x 1 

t = a, 

2 2 



a; 1 
s = — I — a. 
2 2 



Replacing (j2*Uj) in (fT7|) we obtain 



/(x,7i) 



a: — ex — n 



+ n In (x + a) , x > v 2n. 



(20) 



(21) 



We shall now find g(x, n). Using (j2H) in 0, we get 

1 dg dg 1 



2cr (x + a) dn dx (x + a) 



or 



( x _|_ cr) ^ + ® g - 1 
5n cte 2a 



(22) 



Solving (J22J), we obtain 



g(x,nj 



In -2 



x — n + xa 
x 2 — 2n + xcr 



+ C(x + a) 



where C(x) is a function to be determined. Imposing the condition we 
have 

= c/(x, 0) = - ln(-2) + C(2x). 



Thus, 

#(x,7i) 

We summarize our results in the following theorem. 



I In 
2 



1 f- + i 

2 V<r 



(23) 



Theorem 1 In £/ie region x > \j2n. the Hermite polynomials admit the 
asymptotic representation 



H n (x) ~ $i(x,n) =exp 



r i 2 -ffx-n 



+ n In (a + x) 



n — > oo, 



(24) 



where a(x, n) was defined in HQ) . 
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Using the reflection formula (J2J) we can extend our result to the region 
— x > \/2n and obtain: 

Corollary 2 In the region x < —\/2n, the Hermite polynomials admit the 
asymptotic representation 

H n (x) ~ § 2 {x,n) =(-l) n exp 




To illustrate the accuracy of our results, in Figure [T] we graph H^(x) and 
the asymptotic approximations 4) and $2(2;, 4). 



x + ax 



n 



+ n In (cr — x) 



(25) 
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2.2 The transition layer 

We shall now find an asymptotic approximation for |x| ~ y/2n. We will 
consider the ~ \[2/n and find the corresponding result for x ~ —\[2/n 

by using (J2J). From (J23jl we have 



$i(x, n) ~ exp 
We define the function G n (x) by 
if n (x) = exp 



n 3 , — 

— In (2n) rz + y2nx 

2 y ' 2 



x 



2n 



n 3 > — 

— In (2n) n + v 2nx 

2 K J 2 



Using (pj)) in ((TJ) we get 

"n + 1 



In (2n + 2) - | In (2n) - - + ^2 (n+1) - \Phi 



exp 



+V2nG n + G; = 2xG n (x). 
We introduce the stretch variable (3 > defined by 



x 



(26) 



G n +i 
(27) 



x = 

and the function A (/3) defined by 

G n (x) = A 

From (128) ) we have 
"n + 1 



2n+ 4" 

716 



x 



2n) ns 



exp 



In (2n + 2) - | In (2n) - - + ^2 (n + 1) - v 7 ^) 
2n + /5n~e , n — > oo. 



x 



(28) 



(29) 



(30) 



Using (2EJ) in J2HJ) we obtain 



G n+ i(x) = A 



( v 7 ^- V2n + 1 + -4- J + 



(31) 



A(/3)-^A'(/3)n-UiA"(/3)n-§, n - oo, 
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and 



2nG n + G' n - 2xG n {x) 



2nA ((3) + A' (0) - 2(3A (/?) tT* . (32) 



Using (|3*0J). (}3*T]) and (J32*j) in (|2*7jl we obtain, to leading order, the Airy 
equation 

A" (/?) = 2V/2/9A (0) . (33) 



Thus, 



A (/?) = d Ai (v 7 ^) + C 2 Bi , 



(34) 



where Ai (•) and Bi (•) denote the Airy functions and C±, C2 are to be deter- 
mined. Replacing (|2*Hj) and ()34|) in ()26|) we have 



H n (x) ~ exp - In (2ne) + v^/fas Ci Ai (V2/3J + C 2 Bi (V2/3) . (35) 
To find Ci, C 2 we shall match (jHU) with (J2H). Using (|2B jl in pi jl we get 



$i(x, n) ~ exp 



7 

-ln(2ne) + >/2/?nS 

2 3 



_5 1 1 

2 4 n e, 



(36) 



as P — > 0. Using (|28j) and the well known asymptotic expansions of the Airy 
functions 



Ai (x 



2^ 



1 / 2 3 

— — =exp — £2 |i 4 ; x — ► 00 



r>. / \ 1 (2 3\ 1 

Bi (x) ~ — p= exp -x 2 ] x 4 , x — > 00, 

'7T V 3 



in (|33j) we have 
n 



exp 



In (2n) 



-n + v 2nx 



exp 1 1 [1 + In (2n)] + V^/fos } (37) 



and 



d Ai (V2p\ + C 2 Bi 



2*-<A C 2 /2J . 



v 72l/3i eXP V 3 ^ 2 J + 0F2i/? 



(38) 
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Matching to (|37j) and (JHSJ) , we conclude that 

Ci = V27m», C 2 = 0. 



(39) 



This completes the analysis. Combining the results above, we have the fol- 
lowing result: 



Theorem 3 For x ~ \2n, the Hermite polynomials have the asymptotic 
representation 



H n (x) ~ ® 3 (x,n) = exp 



n 3 > — 

— In (2n) n + y2nx 

2 y 1 2 



(40) 



x v27rn6 Ai 



\[2 (x — V2n) ne 



n — ► oo. 



Use of the reflection formula (J2J) provides the corresponding result for 
x rs —y2n. 



Corollary 4 For x ~ — v 2n, t/ie Hermite polynomials have the asymptotic 
representation 



H n (x) ~ $ 4 (a:,n) 



x V27rn6 Ai 



-1) exp 

71 (i 



71 3 / 

— In (2n) 7i — y2nx 

2 K ' 2 



(41) 



71 — > OO. 



2.3 The oscillatory region 

2 

We now study the region bounded by the curve n — y, where the zeros of 
H n (x) are located. In this region, the solution is a linear combination of (|24j) 
and 

H n (x) ~ $5(2;, 71) = n) + K^dx, n), n — > 00 



with |x| < v2ti and K 1 ,K 2 are constants to be determined. We shall require 
$ 5 (x, 71) to match $3(2;, 71) asymptotically in the local variable j3, i.e., it must 
satisfy the limiting condition 



lim$ 5 (/5, n) = lim $ 3 (/?,n). 

/3— >0 /3^— 00 
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Writing (|4U|) in terms of /3, we have 

$ 3 (/3,n)=exp - In (2rae) + -\/2/?ra5 Ai (^2/3 

Using the asymptotic formula 

1 



Ai (x) 
in (02) we get 



sm 



2 3 7T 



-re) - 



-oo 



$ 3 (/3, n) ~ exp - In (2ne) + 



x sm 



3 I 



-oo, 



which can be rewritten as 



$ 3 (/3,n)~2 8n6/3 4 exp — In (2ne) + V2(3m 



x 



| 7 3 \ / 1 7 3 

exp ( --24/32 1 + iexp I 



/3 



-oo. 



- In (2ne) + V2(3n^ - -2^fi 
2 3 



Using in 1(2^1 . we have 
$i(/3, n) ~ exp 
Similarly, using (J2~%|) in (J23j) . we obtain 

$2(/3, n) ~ exp 



5 1 „ 1 



- In (2ne) + V2(3n^ + -25/35 
2 3 



5 1 ^ 1 



where we have used 



1 \ n 

-I) exp 



x + ax — n 



+ nln (<r — x) 



exp 



x + ax — n 



(42) 



(43) 



2-sna/3-i, /3->0. (44) 



2-8n6/3"3i, /3^0, (45) 



Comparing (|43|) with (J44"|) and (J45j) we conclude that K\ — 1 
therefore 

$5(x, n) = $l(x, 77.) + $2(^, 



+ nln (x — a) 

K 2 and 
(46) 
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Since —y2n<x< y2n, we set 

x = v / 2n sin (9) , 
Using (|37)l in (JTHJ), we have 

a = V2n cos (9) i 
Replacing (gBJ) in we get 



7T 7T 

- < 6> < -. 
2 2 



exp 



ar — crar — n 



+ nln (cr + x) 



exp <| — [In (2n) — cos (29)] — n 



- sin (29) + 



and 



2 V a 



ex p (-10 

2 cos (0) ' 



Similarly, replacing (|48j) in (J25j) . we obtain 



;-l) n exp 



x + ax — n 



n 



exp <J - [In (2n) - cos (29)} + n 



+ n In (a — x) 
1 sin (29) + 9 



and 



1 / _x 

2 V a 



2 

ex P (10 

y/2 COS (9) 



Using (gnj)-(j32J in (HSJ), we have 



2nsin (#) , n 
1 



cos (9) 



exp || [In (2n) - cos (20)]} 



x cos < n 



7T 



2 sin (29) + o 



7T 7T 

— < 6 1 < -. 
2 2 



Thus, we have proved the following: 
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Theorem 5 In the region \x\ < v2n, the Hermite polynomials have the 
asymptotic representation 



2nsin (9) 



cos (0) 



exp{|[ln(2n) -cos (20)]} (54) 



x cos < n 



I s in(20) + 0-| 



n — > oo, 



with -f < 6 < f . 



In Figure El we graph 



2n sin (9) exp |-| [In (2n) - cos (20)]} 



and 



cos (0) 



cos < n 



± sin (20) + - | 



+ 



with n = 20. We only include the range < 9 < f , since both functions are 



even. 



The same results obtained in this section were derived in [§| using a dif- 
ferent method, based on the limit relation between the Charlier and Hermite 
polynomials [TH] . 



3 Zeros 

We shall now find asymptotic formulas for the zeros of the Hermite poly- 
nomials using the results from the previous section. Let's denote by > 
Q > ■ ■ ■ > Q the zeros of H n (x) , enumerated in decreasing order. Then, it 
follows from (p>f|) that 



r" 



2nsin(7£) , 



n 



oo 



(55) 



where tJ? is a solution of the equation 



n 



1 sin (2r£) + - \ 



TV 



+ = (1 - 2Jfe) ~, l<A;<n. (56) 
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Figure 2: A comparison of the exact (solid curve) and asymptotic (ooo) 
values of H 2 q (x) in the oscillatory region. 
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Solving (J56|) numerically and using (J55)) we get very good approximations of 
Q . One could also solve (JoT)|) exactly (as we did in 0) and obtain a Kapteyn 
series expansion for t£ 

where iV = 2n + 1 and Jj (•) denotes the Bessel function of the first kind. 
However, (JSTj) is difficult to analyze asymptotically. Hence, we will take a 
different approach and find an approximation for t£ from (j^Sj) trough per- 
turbation techniques. 

We will consider two cases: k = 0(1) which corresponds to the largest 
zeros of H n (x) and k = O f~), related to the zeros close to x = 0. 



3.1 Case I: k = 0(1] 
Replacing 



7T 



T, 



2 

in ([55)1 we obtain, as n — > 00 



^ai(fc)n-t (58) 



1 1 1 _1 K K 3 / 2 \ 

ai = — ac 3 , ao = — K 3 , as = , gu = (k —5 

2 ' 2 ' 3 120 30 V ; 

K~i 43 

a 5 = (3k 4 + 350k 2 + 1400) , a 6 = «. 

5 8400 V ; ' 6 16800 ' 

a 7 = ^j^ty + 350k 4 - 980k 2 - 11200), 

a 8 = --^^(13k 6 + 475k 4 + 1400k 2 + 17500), (59) 
8 63000 v ; v ; 

59 43 

da = K H K 

67200 34496000 

_ 17 

a M = — (23817k 8 + 2608760k 6 - 4592280k 4 

1397088000 v 

- 51744000k 2 - 664048000), 



3 



with 



(/c) = 3tt (4Jfe — 1) . (60) 
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Using (|5BJ)-(|5nj) in (J5SJ), we get 



— i i /■. ■■: _i 1 _1 K 2 + 80 _5 

Cr ~ v2 na n e + -n 2 s-n b 

1 ° 4 6 40«l 



11k 2 + 3920 _3 5k 4 + 96k 2 + 640 _^ 

n 2 -| . n e (61) 



ii 



-n e 



179200 7680k§ 
823k 6 + 647200k 4 - 2464000k 2 - 25088000 

258048000KT 

3064 + 33k 2 _ b \ 

H ra 2 n —>■ oo. 

716800 / ' 

3.2 Case II: fc = O (f ) 

We now set 

I Tl I Tl 

k= L2J +1_J= 2 _a + 1 -^ ^ 

where a = frac (|) (the fractional part of |) and j = 0, 1, 2, . . . . Using (jo^j) 
and 

tf = $>0>-' ( 63 ) 

t>i 

in (jo^j) we obtain, as n — > 00 

*i = & 6a = -|, &3 = I (3 + 16e 2 ) , 

&4 = ~192 ( 3 + 64 ^)' ^5 = ^(15 + 800^ + 10240, (64) 
6 6 = _-J_(l5 + 1600£ 2 + 7424£ 4 ), 

with 

£0') = £(2i + 2a-l). (65) 

Using in we obtain 

>n -i 1-3 3 + 8£ 2 _5 3 + 80£ 2 _z\ 

4 ~ V2£ ( n 2 - -n 2 + -^^n 2 >\ n ^oo. (66) 

In Tablel we compare the exact value of the positive zeros of if 20 (%) with 
the approximations given by solving (JHH|) numerically and formulas (|61j) and 
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Table 1: A comparison of the exact and approximate values for the positive 
zeros of i/ 20 (x) . 



/-n 




mi 


fro 


.24534 


.24536 


.24536 




.73747 


.73751 


.73750 




1.2341 


1.2342 


1.2340 




1.7385 


1.7387 


1.7376 




2.2550 


2.2552 


2.2512 


2.2592 


2.7888 


2.7892 


2.7779 


2.7912 


3.3479 


3.3486 




3.3492 


3.9448 


3.9456 




3.9460 


4.6037 


4.6056 




4.6055 


5.3875 


5.3939 




5.3937 



(}66j) . Note that the biggest error corresponds to the larger zero, where the 
asymptotic approximation (|54jl almost breaks down. 

We summarize the results of this section in the following theorem. 

Theorem 6 Letting > Q > ■ • ■ > Q be the zeros of H n (x) , enumerated 
in decreasing order, we have: 

1. 

r( i «i _i 1 _i k 2 + 80 A 

C? ~ V2 - — n H-n 5 - 1 n 5 I , n -> oo, 

V ° 4 640^3 y 

where k = 0(1) and k (k) was defined in \6U\) . 

2. 



n — > oo, 



/— / i 1 _i 3 + 8£ 2 _5 3 + 80£ 2 _ 7 \ 
C ~ v 2£ n 2 n 2 -| n 2 n 2 

where k = | — a + 1 — j, a = frac (|) and £ (j) was defined in Kffty . 
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